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Abstract 


Airborne observations of the 16-39 y spectra of ten 

oxygen-rich stars with excess emission in the infrared have 

been obtained. The stars show excess emission attributed to 

circumstellar dust grains in the 16-39 y region in the form of 

a broad hump peaking near 18 y and falling smoothly to longer 

wavelengths. The emission is similar in character to the 

emission from the Trapezium region of the Orion nebula indicating 

the grain materials are quite similar in these objects. The 

existence of a feature in the 20 y region is consistent with 

the O-Si-O bending resonance expected for silicate material. 

The lack of any sharp structure in the spectra indicates the 

silicate is in an amorphous, disordered form. A simple model 

of small grains of carbonaceous chondrite silicate material in 

a diffuse circumstellar envelope is shown to give a good 

qualitative fit to the observed 8-39 y circumstellar spectra. 

Comparison of the observed spectra with the model spectra 

2 

indicates the grain emissivity falls as 1/A from 20 y to 40 y. 



I. Introduction 


At wavelengths longward of 8 p, the spectra of many 
cool stars show an excess in flux which has been attributed 
to circumstellar dust emission (Gillett, Low, and Stein 
1968; Woolf and Ney 1969) . The circumstellar nature of 
this excess emission has been confirmed by spatially 
resolving the shells from the star (Zapalla et al 1974, 

McCarthy et al 1977, Sutton et al 1977) and its identi- 
fication with dust supported by the lack of any fine 
structure in high resolution spectra (Gammon et al 1972, 

Treffers and Cohen 1974). Less certain are the chemical 
and physical composition of the dust grains, their temper- 
atures and spatial distribution about the star and the 
wavelength dependence of the opacity of the grains. These 
characteristics are important to an understanding of the process 
of grain formation and mass loss from the stars. Further, 
infrared spectroscopy has shown that the grains around 
oxygen rich stars are quite similar to an important com- 
ponent of the interstellar dust grains. Thus information 
on the composition and opacity of the circumstellar dust 
grains has important applications to the understanding of 
other galactic infrared sources. 

The present study of the 16-39 p spectra of oxygen 
rich stars provides new information about the dust grains 
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in the circumstellar envelopes. If the 10 y emission 
feature in circumstellar spectra is due to an Si-0 stretching 
resonance in silicate grains, as suggested by Woolf and 
Ney (1969) , a second O-Si-O bending resonance is expected 

* 

in the 20 y region. The shape of the 20 y resonance depends 
on the chemical/physical composition of the silicate 
material. For instance well ordered silicates, as are 
found on the earth and Moon, all show fine structure in 
the 20 y region. The instrumentation and observations are 
presented in sections II and III and discussed in section IV. 

A simple shell model of silicate grains is compared 

with the observed spectra in section V. 

II . Instrumentation 

The observations were made with helium cooled Ebert- 
Fastie spectrometers using telescopes aboard NASA aircraft 
operating between 41,000 and 45,000 ft. altitude. The 
two channel spectrometer covers the wavelength range 
16-39 y (Forrest, Houck, and Reed 1976) . An As:Si detector 
was used for the 16-23 y band with a resolution 0.5 y FWHM. Prior 
to 1977 the 21-39 y band used a Ge:‘Ga photoconductor built by 
W. J. Moore with a resolution of 1.2 y FWHM. In February 
1977 this was replaced by a more sensitive Ge:Zn detector con- 
figured to give a resolution of 1.0 y FWHM- The 10 channel spectro- 
meter covers the wavelength range 16-30 y and employs 
10 Si:Sb photoconductors with an Intel 8080 microprocessor 
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based data system which demodulates the signals of 10 
detectors simultaneously. This spectrometer and data 
system are described more fully in a separate paper 
(McCarthy, Forrest, and Houck 1978) . The spectral 
resolution was 0.2 y FWHM for all the observations 
except for the spectrum of VY CMa for which 
AA = 0.5 y FWHM. The spectrometers were used 
in conjunction with the 91 cm telescope of the NASA 
Gerard P. Kuiper Airborne Observatory (KAO) and the 30 cm 
telescope of the NASA Lear Jet (VY CMa only) . The beam 
size of approximately 30" on the KAO and 2.7 T -3.2' on the 
Lear Jet was large enough to include most of the flux from 
these compact sources. Standard sky chopping using the 
oscillating secondary was employed. 

III. Observations 

The program stars are listed in Table 1 along with 
their two-micron sky survey number (IRC, Neugebauer and Leighton 
1969) , their AFGL sky survey number (AFGL, Price and 
Walker 1976) , the spectral and variable type, and an estimate 
of the fractional amount of the total emission which is in 
excess of that expected from the star alone ("bolometric 
excess"). The stars were selected on the basis of their 
brightness at 20 y, the presence of the 10 y "silicate" feature, 
and spectral evidence in the visible of an oxygen-rich atmosphere 
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(0/C > 1) . With the exception of IRC + 10420 
all the stars are luminous, red giant stars 

probably in an advanced stage of stellar evolution. IRC 
+ 10420 is an unusual and luminous GO I star with large 
excess emission from dust (Humphreys et al 1973) . 

OH 26.5 + 0.6 shows evidence of having the thickest 
dust shell of the stars in this group, with strong 
absorption at 10 y and a very cool energy distribution. 

It is thought to be a Mira variable red giant which is 
losing mass rapidly (Forrest et al 1978) . The stars 
are arranged here in order of increasing optical depth 
of the circumstellar material, as judged their fractional 
bolometric excesses (column (6)). The methods used in deriving 
this "bolometric excess" is described more fully in section IV. 

The observations reported here were obtained between January 
1976 and May 1978. The line of sight water vapor column density 

was monitored aboard the KAO and ranged from 5-20 y 

precipitable water vapor. The absence of water vapor 

features in the derived spectra shows that the absorption features 

have been adequately removed by the normalization procedure. 
Shortward of 17 y, absorption due to terrestrial CO 2 
affects the spectrum, and it was necessary to make small 
differential corrections to account for this. For spectral 
calibration, the Moon was observed with the same instruments. 

The Moon was assumed to radiate as a blackbody at a temperature 
appropriate for the phase of and position of the beam on 
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the Moon as judged by the 10 y measurements of Geoffrion 
et al (1960) . For the absolute flux calibration on the 
KAO, ground-based measurement of the 20 y flux from non- 
variable sources were used as follows: May 1976, NGC 7027 

^ F 20y = x 16 W/cm 2 y, Becklin et al 1973); November 1976, a 

-15 2 

0ri(F2Q 1J = 1.43 x 10 W/cm y, Morrison and Simon 1973); June 

—16 2 

1977, y Cep (F 2 q u = 5.8 x 10 W/cm y, Morrison and, Simon 1973); 

—16 2 

May 1978, a Her (-F 2Q ^= 3.7 x 10 W/cm y, Morrison and Simon 
1973) . For the Lear Jet and January 1978 KAO observations 
Mars was used as the flux standard at 20 y using the model 
of Wright (1976) . 

The spectra derived under these assumptions are 
shown in Figures la-le. The flux levels have been normal- 
ized for convenient display. In Table 2 the spectra are 
listed. Column (1) gives the number of the spectrum which 
identifies the spectrum in the appropriate figure. Columns 
2 and 3 give the name of the star and the date of the 
observation. Column 4 gives the spectrometer and telescope 
used and column 5, labeled F 2 Q li , 9^ ves the flux level of 
the horizontal tick mark next to the identification numbers 
in the figures; it is also the approximate flux level 
observed at 20 y. In the cases where the infrared 

excess is small (Table 1) , a considerable fraction of the 
flux observed is from the stellar photosphere. 

To estimate the amount expected from the stellar photosphere 
alone, photometry from .84 y to 5 y was fit with a blackbody 
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of temperature appropriate to the spectral type and 
extrapolated to 20 y. This flux will be denoted F 2Q (star); 
ratio F 20 ^ (star) /F 20 ^ is given in column 6. As an example of 
the probable contribution of the stellar photosphere, the 
blackbodies for RX Boo (spectrum #3; Fig. lb) and y Cep 
(spectrum #7, Fig. la) are included. 

Also included in the figures are ground-based, broad- 
band measurements of the non-variable sources made by- 
other observers. The agreement is quite good in all cases 
except for the 33 y observation of y Cep (Fig. la, spectrum 
#7) reported by Hagen et al (1975) . Their flux level 
is approximately a factor of two higher than was observed 
here, though the 25 y flux levels agree within the errors. 

A second spectrum of y Cep taken in January 1978 from 16-30 y 
(spectrum #8 in Fig. la) agrees quite well in shape and flux 
with the 16-39 y spectrum taken in June 1977 in the region 
overlap. Since y Cep is not variable in the infrared at 
20 y (Morrison and Simon 1973) and has not displayed any 
large variations in visual flux from 1974 to 1978 (Mattei, 
1978) , the discrepancy is probably not due to source 
variability. It is felt that the measurements reported 
here are sufficiently reliable to rule out the factor of 
two excess for y Cep at 33 y reported by Hagen et al (1975) . 
This point could be settled by an independent measurement 
of the 33 y flux level of y Cep though this would be 
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difficult to do from the ground because of the large and 
time-variable opacity of terrestrial water vapor at these 
wavelengths . 


IV. Discussion 

All the stars in this study radiate in excess of what 
would be expected from the star alone in the 16-40 y band 
(Table 2) . The probable source of this emission is 
circumstellar dust heated by stellar radiation. The 
character of this radiation is a smooth, blackbody-like 
continuum at long wavelengths with leveling off or 
downturn shortward of 18 y, most evident in the cases 
where the dust emission dominates (e.g.y Cep Fig. la, 

PZ Cas Fig. 1c, VY CMa and IRC + 10420 Fig. le, NML Tau 
Fig. Id) . As this downturn is more rapid than can result 
from blackbody emission it indicates a resonance in the 
emissivity of the dust grains in this region. The spectrum 
of OH 26.5 + 0.6 (Fig. Id) shows an absorption feature at 
around 18 y which is due to radiation transfer effects in the 
extremely thick shell surrounding the star as discussed 
by Forrest et al (1978) . When the optical depth is large, 
a local minimum in flux is observed where the grain 
emissivity peaks due to the decrease in dust temperature 
with distance from the star. 

The spectrum of RX Boo (Fig. lb) , an M8e semi-regular 
variable is somewhat different than the other stars in that 
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the emission is almost flat out to 20 y and then decreases. 
Taking into account the stellar contribution at these 
wavelengths (Fig. lb, smooth curve) , this indicates the 
dust emission peaks at approximately 19 y as compared to 
the peak at approximately 17.5 y seen in the other stars. 

This could be due to a different composition of dust grains 
around this star. Interestingly, the 10 y spectrum of 
this star (obtained by Gillett and Forrest 1972, not shown 
here) is somewhat atypical of M-type stars in that the excess 
peaks at longward of 9.7 y and decreases less rapidly beyond the 
peak. R Leo (Forrest et al 1975) shows a similar behavior. 

Elitzur (1978) has proposed a mechanism for the 
excitation of the main line maser emission of OH at 1665 and 
1667 MHz which is seen from many cool oxygen- rich red giant 
stars. In the present study, the cool stars with thinner 
circumstellar shells, such as RX Boo, X Her, and R Cas, 
are typical of this class of object. In his model, the 
ground state inversion results from a non-equilibrium 
population of upper rotational levels due to pumping by 
a radiation field which falls more steeply than 1^ « 1/A 5 . 

These upper levels are populated by photons at wavelengths 
near X = 35 y, 53 y, 79 y or 119 y. Of the stars that 
have been studied, the emergent spectra of the optically 
thinnest sources with the hottest apparent dust temperatures 
(i.e. with a small infrared excess given in Table 1) goes 
as Fj ^ (1/A) 4 from 30-40 y. In those stars with thick shells 

the observed flux falls off more slowly than this. Thus it is 



not obvious that radiation excitation involving 35 y photons 

as proposed by Elitzur (1978) will work. Closer to the star 

or at longer wavelengths the dust spectrum will probably be 

steeper but this can be offset by a greater contribution from 

the star whose spectrum is probably a blackbody in the Rayleigh- 
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Jeans limit and therefore goes as (1/X) . 

The emission from the dust alone around the stars 
y Cep, PZ cas, VY Ola and IRC + 10420 from 8-38 y is shown 
in Fig. 2. These stars were chosen because of the large 
contrast between the dust emission and the estimated stellar 
emission at these wavelengths . 

For y Cep, the stellar spectrum and excess emission 
from 8-13.5 y was taken from Russell et al (1975). Those 
authors discussed the possibility that some of the emission 
shortward of 8 y could be from the circumstellar envelope. 

The stellar spectrum in the 16-40 y region was assumed to be the 
continuation of the blackbody fit to the 5.5-8 u data and is 
shown as the smooth curve in Fig. la. The difference between 
the observed flux and this blackbody is attributed to circum- 
stellar emission and is plotted in Fig. 2. Integrated over 
wavelength, the circumstellar emission comprises less than 
.05 of the total emission from y Cep (bolometric excess 
< .05 , Table 1), indicating a relatively thin circumstellar 
shell. The subtractxon process is most uncertain near 8 y 
where the stellar continuum can contribute a major fraction of 
the observed flux. 

PZ Cas was observed from 16-39 y at the suggestion of 
F. C. Gillett (1976) who had obtained an 8-13 y spectrum in 
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November, 1975 showing a very large dust emission feature 
peaking at 9.7 u . Using the broad band measurements from the 
IRC (Neugebauer and Leighton 1969) and AFGL (Price and 
Walker 1976) in conjunction with the 8-13 u and 16-39 y 
spectra, it was found that the energy distribution could be 
fit by a 2500° K blackbody shortward of 5 y plus an infrared 
excess longward of 8 p. The blackbody was assumed to 
represent the stellar continuum and was subtracted to give 
the infrared excess shown in Fig. 2. Integrated over 
wavelength, the infrared excess comprised approximately 
0.16 of the total flux observed (bolometric excess = 0.16, 
Table 1), indicating a moderately thick circumstellar shell. 
Again this process is most uncertain near 8 y where the 
stellar continuum can supply a major fraction of the 
observed flux. Beyond 10 y the stellar continuum represents 
less than 10% of the observed flux. There is an additional 
uncertainty in the ratio of 10 to 20 y emission because the 
photometry and spectrometry were done at different times 
and this source appears to vary in the infrared. (Table 2, 
Price and Walker (1976) ) . 

For the stars VY CMa and IRC + 10420, optical depth of 
the circumstellar shell is large enough to modify the stellar 
spectrum and a different procedure to estimate the infrared 
excess was necessdry. In these cases the total observed 
emission was integrated over all wavelengths . In the absence 
of a circumstellar shell, the spectrum of the star was 
assumed to be a blackbody of temperature appropriate for 
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the spectral type (6000°K for IRC + 10420, 2500°K for 
VY CMa) with the same integrated luminosity as the observed 
luminosity. Any emission above this flux level is certainly 
circumstellar excess. For VY CMa the infrared excess integrated 
over wavelength comprises at least 0.5 of the total flux 
(bolometric excess > 0.5, Table 11 while for IRC + 10420 the 
bolometric excess was at least 0.88. For IRC + 10420, the star 
can contribute at most 2% of the flux observed at 8 y and beyond, 
so no subtraction was necessary and the total observed flux is 
plotted in Fig . 2 . For VY CMa, the star may contribute as 
much as 12% to the flux observed at 8 y but this fraction 
decreases to less than 4% at 10 y and beyond. Since the 
contribution of the star to the observed spectrum is 
uncertain when the optical depth and geometry of the 
circumstellar shell are uncertain, the total observed 
flux is plotted in Fig. 2. It is believed to be primarily 
due to circumstellar emission but there may be some stellar 
contribution near 8 y. 

From the derived spectra of 8-39 y circumstellar emission 
shown in Fig. 2, the similarity of the circumstellar 
dust is apparent. There is a peak of emission near 18 y with a 

smooth continuum falling towards longer wavelengths indicating 
a resonance in the grain material near 18 y. At shorter 
wavelengths there is the familiar and more prominent 9.7 y 
emission feature. The existence of these two resonances is 
in qualitative agreement with the silicate hypothesis which 
predicts a feature in the 10 y region due to the Si-0 
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stretching resonance and a feature in the 20 y region due 
to the O-Si-O bending resonance. 

The smoothness of the spectra indicates the grain 
material is different than well-ordered lunar and 
terrestrial silicates which typically show sharp structure, 
especially in the 20 y region {cf Knacke and Thomson 1973, 
Zaikowski et al 1975) . The shape of these spectra 
is quite similar to that observed in the Trapezium region 
of the Orion nebula (Forrest and Soifer (1976) , Forrest, 

Houck and Reed (1976)), but with a higher apparent temperature. 
This is consistent with the idea that the cool oxygen-rich 
red giant stars are the source of an important fraction of the 
interstellar dust seen in Orion. 

V. Shell Model 
A. The Model 

In order to make a more detailed comparison of the 
observed circumstellar spectra in Figure 2 and the theore- 
tical emission some simple models have been considered. In 
the case where the shell optical depth is not large, the 
theoretical emission is particularly simple because each 
grain is observed without attenuation and the radiation 
field which determines the grain temperature falls as 
1/R beyond a few stellar radii. As there are no sharp 
resonances in the observed spectra it is not possible at 
this time to unambiguously identify the grain material 
which is responsible for the emission. However, the 
existence of the two features near 10 y and 20 y is 
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consistent with silicate emission and this will be used 
as a working hypothesis for the present models. It will 
be shown that using measured grain emissivities from a 
likely silicate material and a realistic density distri- 
bution in the circumstellar envelope gives a good fit to 
the observed optically thin spectra. 

Of the laboratory measurements available at present, 
the silicate materials which appear most similar to the 
observed 8-40 y spectra are the carbonaceous chondrite 

meteorites Vigarano and Murchison measured^ byJPenman (1976) 

(cf . Forrest et al 1976) . This is physically reasonable 

because this type of meteorite is believed to represent 
relatively unprocessed material from the primitive solar nebula. 
Other carbonaceous chondrite materials (Zaikowski and Knacke, 

1975) and an artificially produced amorphous silicate material (Day 

1974) show a close similarity to the 10 y feature but the 
second peak occurs at longer wavelengths than the feature 
observed here and by Forrest et al (1976) . Penman measured 
the reflectance from polished samples and applied the 
Kramers-Kronig principle to derive the optical constants 
from 5 y to 40 y. The emissivities of small spherical 
grains were calculated using Mie theory. For the present 
comparison, the emissivities of equal parts of the Vigarno and 
Murchison material has been taken from Penman (1976) , the 
resultant mass opacity, k, , is shown at the top of Figure 3. 

A 

2 

The peak mass opacity of approximately 2800 cm /gm occurs 

2 

at 10 y and a second peak with approximately 1650 cm /gm 
opacity occurs near 18.5 y. 
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The flux from a grain of mass m^, mass opacity ic^, 
at a temperature T^ and a distance D away from the Earth will be 


m i K X 


B X CT i> 


(1) 


2 

where (T) is the Planck blackbody function (W/cm y 
ster) . If the circumstellar shell is optically thin at 
infrared wavelengths, each grain will be seen independently 
and the total flux observed will be: 


F 


X 



Z m ± B x CT.) 


C2) 


Thus the emergent spectrum is determined by the density of 
grains in the shell versus distance from the star (R) 
and the grain temperature as a function of distance from 
the star T^ (R) . The temperature of a grain will be 
determined by energy balance, i.e. 


E in 


= E 


out 


(3) 


In the present case will be dominated by thermal 

radiation so 


E„„. = 4irm. 
out i 


B x (T t ) ax 


( 4 ) 


If the circumstellar shell is optically thin in the near 
infrared (XXI - 5 y) , E . will be dominated by the stellar 
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radiation field, which falls as 1/R more than a few stellar 
radii from the photosphere, i.e. 

2 


E in “ E o <V R > 


(5) 


where E q is the power input at a distance R - R Q from 
the star. 


Equations (3) - (5) thus define the run of temperature 
versus distance from the star. In order to evaluate 
equation C4) it is necessary to know the grain opacity 
at all wavelengths. From 5-40 y, the opacities from 
Penman (1976) shown in Figure 3 were used. Shortward of 
5 ii, it was assumed the optical constants were constant versus 
wavelength and the particles were small so = Kg (-^4 . 

Beyond 40 y, it was assumed that the imaginary part of the 


dielectric constant e 2 fell as 1/X so that k. 


— , 4011 , 

X>40y ^40y ( — — ) 

At the temperatures of most interest in the circumstellar 

shells between 100 and 1000° K, more than half of the 

Planck emission falls between 5 and 40 y so it is believed 

these approximations are not critical to the results. 


Numerical integration of equation (4) showed that « T 

, 6.1 


4.4 


from T = 160 to 1600° K and E . “ T‘ 

out 


out 

below 160° K. 


The density distribution in the circumstellar shell 
is less certain. The mass loss rate from a star is given by 


M = 4ir R z p (R) v (R) 


( 6 ) 


where p (R) is the total density and v (R) is the out 
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flow velocity. Since both the gas and the dust mass loss 

rates have been observed to be fairly constant from stars of 

this type (cf. Forrest et al 1975, Weymann 1963) and the 

circumstellar absorption lines (Deutsch 1960, Weymann 1963) 

and OH emission lines (Wilson and Barret 1972) indicate the 

outflow velocity is constant with time, it is expected that 
2 

p « 1/R in the outer parts of the envelope. Since dust grains are 

not expected to condense where the density is very low and the 

dust is momentum coupled to the gas (Gilman 1972) , one also 

2 

expects the dust density = 1/R in the outer parts of 
the envelope. Close to the point where dust grains 

begin to condense neither the velocitv nor 
rhe fractional amount of dust is well known.' Fortunately, 
it will be shown that at the wavelengths of interest 
beyond 8 u, the shape of the emergent spectrum doesn't 
depend strongly on the density distribution near the region 
of condensation. The point at which dust begins to con- 
dense will be denoted R q X r*, where r* is the stellar 

radius. Theory (Jones and Merrill 
1976)' and observation (Zappala et al 1974) indicate 

r /■£ X 3 for the type of stars considered here. For the purposes 

of this calculation, it has been assumed that all the 

dust condensation occurs at R (cf . Menietti and Fix 1978) . 

o 

Then if radiation pressure on the grains determines the 
flow velocity, the dust velocity will be of the form 

v a (R) = v a (R 0 ) + (1- V /2 


(7) 
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where (R q ) is the initial velocity and v^(R q ) + is 
the final outflow velocity. Two cases have been consi- 
dered: v^ (r ) >> v^, which leads to a a 1/R 2 shell 

structure, and V. (R ) = 0, for which p , a 1 

do a 2 1/2 

R (1-R 0 /R) V 

2 

is infinite at R q but quickly approaches 1/R beyond 

2 R . 
o 

2 

In a preliminary run of the model with p,a 1/R 

a 

and T q = 500° K, an interesting effect was discovered 

when the model was run out to outer limits of 10 R q , 

100 R , and 1000 R . Though 99% of the total flux comes 
o o 

from within 100 R q/ this region contributes only 90% of 
the flux at 40 y and approximately half the flux at 80 y. 
Even though the temperature is low 35-76° K) in the 
outer layers of this shell, the mass is large, and this 
region dominates the emission at longer wavelengths. 

Further from the star than approximately 1000 R q , two 
effects are important in determining the emergent spectrum 
beyond 100 y. First, heating by the interstellar radiation 
field will set a minimum temperature of the grains. With 


the grain emissivity used here, a grain radius of 0.1 y, 

Q a bs = 1' an< ^ t ^ le interstellar radiation field given by 
Allen (1973) , the minimum temperature is approximately 30° K 
in the galactic plane. Second, when a mass of interstellar 
matter equal to that in the circumstellar shell has been 

swept up, the structure of the shell will be significantly 

2 

different than the p a 1/R assumed here. For a space 
density of 1 H atom/ cm , an outflow velocity of 10 km/s. 



and a mass loss rate of 10 ^M Q /yr this happens after 

5 4 

approximately 3 x 10 years at a distance R > 10 R q from 

the star. Since we are interested in the emergent spectrum 

from 5-40 y, these effects will not be important and it 

is sufficient to integrate out to 1000 R q . 

With the density in the shell fixed, the only free 

parameter in the model is T , the dust temperature at the 

distance R q from the star where dust forms. This temperature 

has been varied in order to fit the observed ratio of 

F 10v/ F 18u for t * le s ^ e ^-^ s around PZ Cas and y Cep (Fig. 2) . 

The temperatures are given in Table 3 and the resulting 
spectra are compared to the observed spectra in Fig. 3. 

Also given in Table 3 is the temperature a thin isothermal 
shell would have in order to give the observed ratio 

F 10y /F 18y* 


B. Discussion 


From Fig. 3 it is seen that the models give a good 
fit to the observed spectra. In particular, the shape 
and position of the 10 y and 18 y emission features are 
reproduced. This provides further support to the hypothesis 
(Woolf and Ney, 1969) that silicate grains are responsible 
for much of the excess emission from oxygen-rich red giant 
stars. It also indicates that the grains around these 
stars are probably made of a disordered, amorphous form of 
silicates similar to that found in the carbonaceous 


chondrite meteorites. 
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From Figure 3 it is seen that there is little 
difference between the emergent spectra for the two 
density distributions considered here. This shows that the 
emergent spectrum is not particularly sensitive to 
the density distribution near R = R Q , where it is least 
well known. Thus it is felt that the discrepancies 
between the model and observed fluxes shortward of 9 y 
and longward of 20 y are not due to density effects 
alone. 

Short of 9 p, the models predict more flux than is 
seen. This could be due to either lower emissivity of 
the circumstellar grains at these wavelengths or an over- 
estimate of the stellar continuum which was subtracted to 
derive the excess emission. For both PZ Cas and y Cep, 
the 8 to 10 y spectrum had the approximate shape of the 
model spectrum before subtracting the stellar contribution. 
This indicates that the circumstellar silicates probably 
' have a sharper feature than the silicates measured by 
Penman (.1976) . However, some of the observed 6-8 y 
continuum could be due to the hot grains in the inner 
part of the shell. A measurement of the size of the 
emittxng region at different wavelengths could clarify this 
point. 

Longward of 20 y, the models predict less flux 
than is observed. This indicates that the emissivity of 
the circumstellar grains fall less rapidly than the 
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silicates measured by' Penman C1976) , The observed 
spectra would be fit by an emissivity which fell as 
1/X^ from 20-40 y rather than the l/A^*® dependence of 
the Vigarano and Murchison mixture used here. 

VI. Conclusions 

The study of the 16- 39 y emission from oxygen-rich 
red giant stars and IRC + 10420 with excess emission at 
10 y has revealed: 

(1) The stars show excess emission in the 16-39 y 
region in the form of a broad hump peaking near 18 y 
and falling smoothly to longer wavelengths. 

(2) Except for possibly RX Boo, the emission from 
these stars is quite similar in character and qualitatively 
like the emission seen in the Trapezium region of the 
Orion Nebula (Forrest et al (1976)). This indicates the 
grain materials are quite similar in these objects. 

(3) Lack of further structure in the 20 y region 
prevents an unambiguous identification of the circumstellar 
material. The observed spectra are consistent with 

emission from small grains of silicate material similar 
to the carbonaceous chondrite meteorites in a diffuse 
circumstellar envelope. Comparison of the spectra with a 
simple model indicates that the grain emissivity falls as 
1/X 2 from 20 to 40 y and suggests that some of the observed 
6-8 y continuum could be due to emission from hot grains 
in the inner part of the circumstellar envelope. 
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(4) The spectra of y Cep do not show the 33 y 
emission feature reported by Hagen et al (1975) . 
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Table 1 


(1) 

Name 

(2) 

IRC 

(3) 

AFGL 

(4) 

Spectrum 

(5) 

Variability 

(6) 

Bolametric Excess* 

a Ori 

+10100 

836 

M2 lab 

SRc 

s 

X Her 

+50248 


M6e 

SRb 

s 

y Cep 

+60325 

2802 

M2e la 

SRc 

s 

RX Boo 

+30257 

1706 

M7e-M8e 

SRb 

s 

R Cas 

+50484 

3188 

M6e-M8e 

M 

s 

NML Tau 

+10050 

529 

M6e-Ml0e 

M 

m 

PZ Cas 

+60417 

3138 

M3 la + 

SRa 

m 

VY CMa 

-30087 

1111 

M5e lb pec. 

Lc 

1 

IRC + 10420 

+10420 

2390 

F8-G0 I ++ 


1 

OH 26.5 + 0.6 


2205 



1 


The spectral type and variability type are taken from Kukarkin et at (1969) except as noted: 

+ Humphreys 1970 
++ Humphreys et al, 1973 

*The bolcmetric excess is the fractional amount of the total luminosity which appears in excess 

of that expected from the star alone as described in the text. 

s = small = 0 to 0.05 

m = moderate = 0.05 to 0.5 
1 = large = 0.5 to 1.0 




Table 2 


(1) 

(2) 

(3) 

(4) 


(5) 



(6) + 

Spectrum, 

figure 

Name 

Date 

Note* 

F. 

20u~ 

F 

r 20y 

(star) 







(W/cm y) 

F 

20u 

1, 

la 

a Ori 

Nov. 1976 

2, 

KAO 

1.45 

X 

10" 15 


0.41 

2, 

lb 

RX Boo 

May 1976 

2, 

KAO 

3.8 

X 

10" 16 


0.32 

3, 

lb 

RX Boo 

June 1977 

2, 

KAO 

3.9 

X 

io - 16 


0.31 

4, 

lb 

RX Boo 

May 1978 

10, 

KAO 

4.0 

X 

10" 16 


0.30 

5, 

lc 

X Her 

June 1977 

2, 

KAO 

2.15 

X 

10" 16 


0.33 

6, 

lc 

X Her 

May 1978 

10, 

KAO 

2.1 

X 

10“ 16 


0.33 

7, 

la 

y Cep 

June 1977 

2, 

KAO 

5.8 

X 

10- 16 


0.17 

8, 

la 

y Cep 

Jan. 1978 

10, 

KAO 

6.4 

X 

10" i6 


0.16 

9, 

lc 

PZ Cas 

Nov. 1976 

2, 

KAO 

2.9 

X 

io " 16 

< 

0.1 

10, 

lc 

PZ Cas 

June 1977 

2, 

KAO 

3.7 

X 

io - 16 

< 

0.1 

11, 

Id 

R Cas 

Jan. 1978 

10, 

KAO 

6.8 

X 

io " 16 


0.22 

12, 

Id 

NML Tau 

Jan. 1978 

10, 

KAO 

1.28 

X 

io - 15 

< 

0.1 

13, 

Id 

OH 26.5 + 0.6 

June 1977 

2, 

KAO 

1.08 

X 

io -15 

< 

0.1 

14, 

le 

VY CMa 

Jan. 1976 

2, 

Lear 

6.8 

X 

io - 15 

< 

0.1 

15, 

le 

VY CMa 

Nov. 1976 

2, 

Lear 

7.7 

X 

10“ 15 

< 

0.1 

16, 

le 

VY CMa 

Sept. 1977 

10, 

Lear 

8.0 

X 

io - 15 

< 

0.1 

17, 

le 

IRC + 10420 

May 1976 

2, 

KAO 

2.25 

X 

io -15 

< 

0.1 

18, 

le 

IRC + 10420 

May 1978 

10, 

KAO 

2.0 

X 

io " 15 

< 

0.1 


* Notes (Column 4) : 

2=2 detector 16-39 y spectrometer 
10 = 10 detector 16-30y spectrometer 
KAO = NASA C-141 KAO 91 cm telescope 
Lear = NASA Lear Jet 30 cm telescope 

+ An entry < 0.1 in column 6 means the star contributes less than 10% of the observed flux at 20 y. 
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Table 3 


F 10y / F 18y 

(Observed) Thin Shell 


T 

pd= P ° /r2 P<5= r 2 (1-r°/e, V* 


PZ Cas 
y Cep 


2.95 

6.9 


270°K 

470°K 


400°K 

1250°K 


340°K 

860°K 
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Figure Captions 


Fig. la. The observed 16-39 u spectra of the oxygen rich 
supergiants a Ori and y Cep. The flux calibration and other 
relevant data are given in Table 2. In this and the other 
figures statistical error bars representing ± 1 sigma of 
the mean are shown where large enough to plot. The smooth c ur ve is 
a 4200° K blackbody extrapolated from the 4-8 y data of 
Russell et al (1975) and represents an estimate of the 
probable contribution of the star to the total observed 
flux from y Cep (Spectrum #7) as described in the text. 

The large open symbols are ground-based broad-band 
observations from the following sources: 

squares (Q) Morrison and Simon (1973) 
diamonds (0) Hagen, Simon and Dyck (1975) 
triangle (Z\) Low, Rieke and Armstrong (1973) 

The lower diamond at 33 y refers to the y Cep spectrum below 
it (Spectrum #7) . 

Fig. lb. The observed 16-39y spectra of the oxygen rich 
red giant star RX Boo. The flux calibration and other 
relevant data are given in Table 2. The smooth curve is 
a 2500° K blackbody extrapolated from shorter wavelengths 
and represents our estimate of the stellar contribution 
to the spectrum just above. The open square CO) is the 
ground based 20 y measurement of Morrison and Simon (1973) . 

Fig. 1c. The observed 16-39 y spectra of PZ Cas (M3 la) 
and X Her (M6e) . The flux calibration and other relevant 
data are given in table 2. The open square (□) is a 



broad band 20 y measurement of X Her by Morrison and Simon 
11973). 


Fig. Id. The 16-39 y spectra of the Mira variables R Cas, 
NML Tau (= IK Tau) and OH 26.5 + 0.6. The flux calibration 
and other relevant data are given in Table 2. Some of the 
OH 26.5 + 0.6 data here has been reported previously by 
Forrest et al (1978) . 

Fig. le. The 16-39 y spectra of the super giant stars 
VY CMa and IRC + 10420. The flux calibration and other 
relevant data are given in Table 2. Broad band observa- 
tions by Hagen et al (1975, diamonds) of VY CMa and 
Humphreys et al (1973, upside down triangles) of 
IRC + 10420 are included. 

Fig. 2. The excess emission from y Cep, PZ Cas, VY CMa 
and IRC + 10420 as described in the text. 

Fig. 3. Model fits to the excess emission from PZ Cas and 

y Cep (solid lines) as described in the text. The crosses 

2 

are for a density distribution p^= P Q /R and the diamonds 



for a density pj= 1 . The assumed inner 

a ° “2 172 

R^Cl-R^/R) 1 ^ 

temperatures are, from the top 400° K, 340° K, 1250° K, and 

860° K as noted in table 3. The squares are the opacity 
2 

(k , cm /gm) of the meteorite mixture used in calculating 

A 

2 

the emergent spectra. The peak opacity of 2800 cm /gm occurs 
at 10 y. 
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Figure lb 




( arbitrary normalization) 
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Figure 13 



F ^ (arbitrary normalization ) 
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